WAVELET FRAME BIJECTIVITY ON LEBESGUE AND HARDY SPACES 



H.-Q. BUI AND R. S. LAUGESEN 
Dedicated to our friend and mentor Guido Weiss. 



Abstract. We prove a sufficient condition for frame-type wavelet series in LP, the Hardy space H^, 
and BMO. For example, functions in these spaces are shown to have expansions in terms of the Mexican 
hat wavelet, thus giving a strong answer to an old question of Meyer. 

Bijectivity of the wavelet frame operator acting on Hardy space is established with the help of new 
frequency-domain estimates on the Calderon-Zygmund constants of the frame kernel. 

1. Introduction 

One of the first wavelets ever considered was the Mexican hat (see Figure [I]) , which is 

- (1 - a;2)e-^'/^ xeR. 
Morlet |29j generated the Mexican hat system 

^ i'jAx) = vl'^^Vx ~k), j,ke Z, 



d and applied it to geophysical problems as if it were orthonormal (which it is not). Daubechies |10l 

p. 75] explained why Morlet's approach succeeded: the system is an almost-tight frame, with its frame 
I— I operator 

being a bijection on , and indeed being rather close to the identity. In other words, analysis followed 
2^ by synthesis with the Mexican hat system provides almost-perfect reconstruction in Li^ . Hence the 

frame operator is surjective, and so every square integrable function can be decomposed into a norm 
convergent series in terms of the wavelet system ipj.k- That is, each g has an unconditionally 

convergent series expansion 



T— I j,k 

for some coefficients Cj^k- Do such series expansions hold also when g G ior p 2? 
• '~j Meyer highlighted the extent of our ignorance on this issue when he remarked that "we do not know 

rN whether the functions 2^ /'^ip{2^x-k),j, k eZ, form a complete set in LP {R) fori <p<oo" gTlp. 137]. 

In other words, are finite linear combinations of the tpj^k dense in L^l Note that proving existence of 

series expansions in would be even stronger than completeness. 

Two counterexamples reveal the subtlety of the problem. Tao [3T] showed for each p G (1,2) that 

there exists a smooth, compactly supported synthesizer ip that is "good" on (generating a wavelet 

frame there) but is "bad" on (having non-bijective frame operator). For p £ (2,oo), Tchamitchian 

and Lemarie [27l pp. 130,136] found examples with even worse behavior. Their system {ipj,k} fails even 

to be complete in L^. 

Despite these obstacles, we recently answered Meyer's completeness question in the affirmative, by 
proving L^-completeness criteria for general wavelet systems [2 [5] , and verifying those criteria for the 
specific Mexican hat example. Yet the stronger question of series expansions in remained open. 

This paper answers the stronger question, by proving criteria for wavelet series expansions of arbitrary 



functions in L^, 1 < p < oo, as well as in the Hardy space and BMO (see Theorem 1.1). These 
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bijectivity criteria on the frame operator are explicitly computable. We verify them numerically for the 
Mexican hat system, in Section [6j 

To prove the main results (Theorem |l.l| and Corollary 1.2 ) we develop explicit norm estimates on gen- 
eral Calderon-Zygmund operators (Proposition 3.11 and computable bounds on the Calderon-Zygmund 
constants for the wavelet frame kernel (Theorem 4.6). 

An innovative feature of our work is that we estimate the wavelet frame kernel in the frequency 
domain. Prior approaches proceeded in the time domain: see [211 Ch. 5], [27, p. 165], and related 
approaches [551 [32]. Estimates in the frequency domain are better for our purposes, for three reasons. 
First, we can estimate the Calderon-Zygmund constants while keeping the analyzer and synthesizer 
together instead of splitting them apart (as discussed in Section 4.5). Second, in the frequency domain 
we need not preserve cancelations, because V' is nonnegative for many examples (such as the Mexican 
hat in Figure [I]). Third, the construction of a "good" analyzer is easier in the frequency domain, where 
in order to obtain almost-perfect reconstruction of LP' we need only satisfy a discrete Calderon condition 
and ensure that the "overlaps" are sufficiently small (see Section [6]). 

The techniques in this paper are largely unrelated to our earlier work [Il[l], which relied on L^-type 
methods in the frequency domain and did not use the Calderon-Zygmund theory. 



Wavelet definitions. To state the main result, we need some definitions. Fix A, _B € M with 

\A\ > 1 and B > 0. 
Call A the dilation factor, and B the translation step. We rescale functions t/j,(j) E LP by 

iPj^kix) = \A\^/^i!{A^x - Bk), (f>j^k{x) = \A\^^^<l>{A^x - Bk), j, k e Z. 
The wavelet synthesis operator associated with ^ is the map 

c = {cj,fc} ^ S{c) ^ Cj^k^j.k 

where the coefficients Cj^k are complex numbers. The analysis operator associated with is the map 

f^T{f) = {B{!,<P^^k)L- ),.k& 

where {f,g)L^ — J^fgdx. (The "i?" in the definition of analysis appears only for later convenience. 
Some authors use Slog \A\ instead, to emphasize their Riemann sum interpretation with respect to the 
continuous parameter case; see [lOj.) Sometimes we write and T^, to emphasize the dependence of 
these operators on the underlying synthesizer and analyzer. 

The wavelet frame operator consists of analysis followed by synthesis: 



ST{f)=B ^ (/,0,-fc)L^V',,fc- 
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Main results. Write X{£_) = £_ for the identity function on M, and W''^^ = T4^'='P(M) for the usual 
Sobolev space. Also, write / < g to mean / < (const. )(7. 

Two quantities Mi{iIj, 0) and Al^oiip, 4>) will be defined later, in Sectionjs] Our main result says that 
if these quantities are less than 1, then the wavelet frame operator is bijective on and BMO and 
hence on for 1 < p < oo. 

Theorem 1.1 (Invertibility of frame operator). Assume ip, cj), Xip, X(j> e W^'^ W^'^ with tp{0) = 
(f>{0) = 0, and that the derivatives decay near the origin and infinity according to 

i0'(ouv^'(oi<|j^jj._5/2^ 1^1 ; I; (1) 

for some e > 0. 

If Mi{%lj,(j)) < 1 and Mao(4',4') < 1, then ST extends from to a bounded, linear and bijective 
operator on H^, on for p G (l,oo), and on BMO. 

The proof, in Section [5j involves a quantitative perturbation estimate for wavelet frames: we show ST 
lies within distance 1 of the identity operator on and hence is invertible there. Clearly this approach 
depends sensitively on the choice of norm for . We employ the L^-atomic norm, as defined in the 
following section. Of course, after the frame operator has been proved bijective, it remains bijective 
under any equivalent norm. 

The constants Mi and Moo hi Theorem |1.1| depend on two types of information: (i) Calderon- 
Zygmund bounds on the frame operator, which we control with frequency-domain estimates on the 
analyzer and synthesizer, and (ii) L^-norm bounds on the frame operator, which we obtain in practice 
from Daubechies-type frame conditions. (The original Daubechies condition can be found in .lU] . A 
later variant is due to Casazza and Christensen [6j.) The precise dependence of Mi and Moo on these 
quantities is explained in the definitions in Section[5j and then revisited in a specific situation in Section[6] 
when we treat the Mexican hat. 

Next we give a strong answer to Meyer's completeness question about non-orthogonal wavelets, for 
the class of wavelets in our main theorem. 

Corollary 1.2 (Surjectivity of synthesis). Assume all the hypotheses of Theorem \l . l\ hold. 

Then every function / in , or BMO, or LP for 1 < p < oo, can be expressed as a norm convergent 
wavelet series of the form f — Cj^k4'j.k- 

We apply the corollary in Section [6] to the Mexican hat example. There we construct an analyzer 
which, when combined with the Mexican hat synthesizer, generates a frame operator lying less than 
10"'^ from the identity on and less than about 10~^ from the identity in H^. Thus the frame operator 
is a bijection, and so the Mexican hat wavelet can reproduce each function in H^, and BMO. 

Note that wavelet systems cannot reproduce L-'^-functions with nonzero integral, because the wavelet 
ip and all its dilates have integral zero. Interestingly, reproduction in with < p < 1 can be achieved 
by using only the small scales (terms with j > 0) , as shown in [5S1 §4.4] . 

Open problems that a motivated reader might pursue include higher dimensional analogues of The- 
orem [lT] and the case of spaces for p < 1. 

Relevant literature. Motivating literature for our approach includes: the frame work of Chui and 
Shi [5] , the phi-transform theory of band limited dual frames for Triebel-Lizorkin spaces by Frazier and 
Jawerth [TSlIinillT] with recent extension to irregular translations by Cabrelli, Molter and Romero [S], 
the co-orbit theory of Feichtinger and Grochenig [T31 [TH [53] , and the non-band limited but oversampled 
"approximate duals" of Gilbert, Han, Hogan, Lakey, Weiland and Weiss [18]. Note that Gilbert et al. 
treated translations in a lattice for the full scale of Triebel-Lizorkin spaces. Their work was extended to 
irregular translation grids (still oversampled) by Li and Sun in the restricted case of L^, 1 < p < oo. 

Underlying the approaches of Feichtinger and Grochenig and Gilbert et al. (and Frazier and Jawerth 
[TOl §4], Bui and Paluszyhski |4j, and Li and Sun ^26j) is the fact that a sufficiently dense oversampling 
of translations and dilations must yield approximate reconstruction in any reasonable function space, 
for any reasonable analyzer and synthesizer. The reason is that by the Galderon reproducing formula, 
perfect reconstruction holds in the limit of infinite oversampling (the case of continuous parameter 
translations and dilations). Unfortunately, those papers do not specify the degree of oversampling 
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required to guarantee their conclusions. In other words, the dilation factor A must lie sufficiently close 
to 1 and the translation step B sufhciently close to 0, but how close would constitute "sufficiently close" 
is simply not specified. 

Theorem 1 1 . 1 1 avoids oversampling completely, by constructing a criterion that is usable for fixed values 
of A and B. Further, the theorem relaxes the band limitation of Frazier and Jawerth, so that one may 
treat non-band limited examples such as the Mexican hat. The cost of that relaxation is merely that the 
frame operator yields approximate rather than perfect reconstruction: our frame operator is bijective 
rather than equalling the identity. 

Let us now contrast our work with the established theory of orthonormal wavelets. In that setting, 
the frame operator gives perfect reconstruction on and is bounded on other function spaces such as 
and L^, assuming reasonable smoothness and decay of the wavelet. The frame operator therefore 
equals the identity on those spaces, simply by density of , and thus one obtains wavelet expansions 
in whole families of function spaces [241 127) . In this current paper, we do not assume orthonormality or 
perfect reconstruction on . We assume instead that the difference between the frame operator and 
the identity has small norm in . The challenge in Theorem |l.l| is to show (under suitable hypotheses) 
that the norm of this difference remains small when considered in the Hardy space and BMO, so that 
the frame operator is bijective on those spaces. 

Our work is unrelated also to multiresolution analyses that generate pairs of perfectly reconstructing 
analyzers and synthesizers, such as (7, jT, 30_ . Such methods do not apply to the Mexican hat wavelet, for 
example, because it satisfies no scaling or refinement equation. Theorem |1 . 1 1 takes a different approach: 
it finds sufficient conditions for a given analyzer and synthesizer to generate a frame-like structure. 

2. Definitions 

We employ the Fourier transform with 27r in the exponent, 

so that Parseval's identity says {f,g)L^ = {f,g)L2. 

Next we specify the norms to be used on Hardy space and BMO, and define the Calderon-Zygmund 
operators that will act upon those spaces. As always, we write ~ LP{M.). 



Atomic Hardy space, and BMO. We work with the Hardy space = H^{M.) defined by atoms, 
as follows. Call a g an atom if a is supported in some bounded interval / with ||a||^2 < |/|^^/^ and 
Jja{x)dx = 0. Notice ||a||Li < 1. Define 

oo 

= { / e : / = XjOj for some atoms aj and 



constants Xj € C satisfying ^ |Aj| < oo}. 



3 = 1 



(Note the series / — J^j ^j^j converges in .) Then is a Banach space under the atomic norm 

oo oo 

||/|Ui=inf{^|A,|<cx):/ = ^A,a,}, 
and it imbeds continuously into with 

ll/IU^ < ll/llff- 

The dual of the Hardy space is (H^)* ~ BAIO, the space of functions with bounded mean oscillation. 
We employ the dual space norm on BMO: 



\g\\BMO 



= sup{|5(/)| : / e with ll/llffi = 1}. 
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Calderon Zygmund operators (CZOs). 

Definition. Call Z a generalized Calderon-Zygmund operator if Assumptions la, lb, Ic, 2 and 3 below 
all hold, for some constants Ci, C2, C3 > 0. 

Suppose K{x, y) is a measurable, complex valued function on M x M. 

Assumption la. 

\K(x,y)\ < Ci-. r whenever x^y. (2) 

Assumption lb. 



\x — Xq\ 1 

\K{xo,y)~K{x,y)\<C2 . _ ,3 whenever |a; - xol < -|a; - y|. (3) 

\x y\ Z 



Assumption Ic. 



\Kix,yo)-K{x,y)\<C3\!^^^ whenevev \yo - y\ < -\x - y\. (4) 



\yo-y\ , I 1^1 

\x-y\2 whenever \yo - vl < ^ 

Assumption 2. Z : — > is a bounded linear operator such that if / G has compact support 
then 

iZ.f){x)^ / K{x,y)f(y)dy whenever a; ^ supp(/). (5) 

Assumption 3. If / G has compact support and integral zero, then Zf e L^OL'^ has Jg^{Zf ){x) dx — 0, 
and Z*f€L^nL^ has J^{Z*f ){x) dx = 0. 

Note that the integral representation ([5| need not hold when x S supp(/), for example if Z is the 
identity operator and K = 0. 

Remark. This definition of generalized CZO follows the treatment by Hernandez and Weiss [24]. A more 
general distributional definition was developed earlier by Coifman and Meyer [28] ; see also Grafakos [21] . 

A function used later. Define 



so that D{() - 7/C as C ^ 00. 

3. CZOs on atomic Hardy space 

Our wavelet results depend on boundedness of generalized Calderon-Zygmund operators. Bounded- 
ness of CZOs on H^,LP and BMO is well known [551 HIl IUji but boundedness is not enough for our 
purposes. We need a certain operator norm less than 1, in order to prove invertibility of the wavelet 
frame operator by using a Neumann series. Thus we must have explicitly computable norm bounds on 
CZOs. The point of the next proposition is to prove such explicit bounds. 

The constant used below is 

dp) = (2^./(.+i)fcl))<^-^'^-^^^/^^ . ^ as p M. 
V p — I / P — ^ 

Notice c(2) = 1. Write p' for the conjugate exponent: ^ + ^ = 1. 

Proposition 3.1. Assume Z : ^ LF' is a generalized Calderon-Zygmund operator with constants 
Ci,C2,C3. FixC>2. 

Then Z has a bounded linear extension to , to for each p e (1, 00), and to BMO, with norms: 
\\Z\\m^m < 2v^ \\Z\\l-^^l-^ + DiQCs, (6) 
\\Z\\lp^l. < c{p){^\\Z\\l2^l- + ^^Csf^^^-'WZWl-^^/^^ when l<p<2, (7) 

\\Z\\lp^lp < c{p'){./32C\\Z\\l2^l- + T^C2Y'^^^''^\\zf/.r^^, when 2 < p < 00, (8) 



\\Z\\bmo^bmo < 2^\\Z\\l2^l-+D{C)C2. (9) 
Equality holds for p — 2. 
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The purpose of including the free parameter C. in the proposition is to improve the bounds by allowing 
trade offs between the size of \\Z\\]^2_^i^2 and of C2 and C3. For instance, ^ = 50 is a nearly-optimal 
choice for the Mexican hat example in Section [6j whereas choosing C = 3 would give results an order of 
magnitude worse. 



In order to prove Proposition 3.1 at the end of this section, we will find explicit bounds saying 
(roughly) that generalized CZOs map atoms to molecules, and that molecules decompose into atoms, so 
that molecules belong to the Hardy space. The next two lemmas phrase these known results in a way that 
avoids the molecular norm, so as to prove the proposition more directly and obtain computable constants. 
Note that here, and in later proofs, we must not switch between different types of atoms, because doing 
so could introduce large additional constants into the norm estimates. We must consistently use 
atoms. 

The first lemma says that Z maps each atom to a function that decays like the square of the distance 
from the support of the atom. 

Lemma 3.2 ("CZOs map atoms to molecules"). Let C > 3. If Assumptions la, Ic and 2 hold and a{x) 
is an atom supported in a bounded interval I centered at yo G M, then Za G O L"^ with 

\iza){x)\<-^^, xeiciy, 
[x-yor 

where the constant is 



1 / cnC^ + 3) ^ 

Here (I denotes the interval having the same center as / and C times the length. 

Proof of Lemma \37^ We adapt parts of the standard proof; see for example [211 Theorem 5.6.8]. 
First, Za G by Assumption 2, since the atom a belongs to L^. 
Next, if X lies outside the support / of the atom a then 



{Za){x) = J K{x,y)a{y)dy by Assumption 2 

[K{x,y)-K{x,yo)]a(y)dy (10) 



since jj a{y)dy ~ 0. 

Now suppose ?/ G / and x G [CIY- Then 



|yo-2/|<^|/|<^|/|<^|x-y|. 



where the middle inequality uses that C ^ 3. Assumption Ic then implies with ( 10 1 that 



{x - yY 



\{Za){x)\<C, I ^i^^^^\a{y)\dy (11) 



(„p y)2 ,1/2 

<C3( / -, ^ dy\ ||a||L2 by Cauchy-Schwarz 



by a change of variable, where 



y - 2/0 , x-yo 

z — - — -- — and w = - — — ; — . 
|/|/2 |/|/2 
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Evaluating the last integral shows that 

|(Za)(i)| <C3,/^l/ ';4 ii°iU' 




since the atom a has L^-norm at most |/|^^/^. Note Iwj > C (since x ^ a-nd that the function 
w I— ?> w^{'w'^ + 3)/(w^ — 1)3 is decreasing for w > 1. Hence we may replace w with inside the square 
root in the last formula. The lemma now follows by substituting the definition of w into the remaining 
factor of 

Notice Za is locally integrable because it belongs to L^, and is globally integrable because it decays 
like a;~^ at infinity. Thus Za e L'^. □ 

Lemma 3.3 ("Molecules belong to Hardy space"). Suppose M £ and I is a bounded interval centered 
at yo e M. Let C4 > and C > 3. 
If 

\Mix)\ < xeiCir, (12) 

and Jjg M{x) dx — 0, then M € with norm 

\\M\\hi < 2y/C\T\\\M\\L2 + ^(%/57+ v/320/3)C4. 

The proof adapts the essential ideas from Coifman and Weiss's decomposition of a molecule into a 
sum of atoms Theorem C]. They used dyadic scaling, 7 = 2, whereas we obtain a better constant by 
optimizing over 7 > 1 in the proof below. Also, they chose the initial radius R to depend on the L"^ 
norm of M, whereas our choice of R will depend on C and /. But the central idea of decomposing M 
into atoms supported on annuli comes directly from their work. 



Proof of Lemma \3.,% We may suppose M ^ 0, and that the interval is centered at yo = 0, by a trans- 
lation. Let 7 > 1. (We will later choose 7 — 4.) Define 

._C|/| 

= h7^i?,7^i?], J =0,1,2,3,... 
f}o = /o = [-i?,i?]-a, 

% = I, \ = [-7J"i?, 7^"i?] \ [-7^-^i?, 7^"'i?], j = 1, 2, 3, . . . 
Ij In, , 

Clearly Mj is supported in C Ij, and Mj is square integrable and has mean value zero. 
We break the proof into several steps. 

Step 1. We will prove Mj /Xj is an atom, for numbers Xj that satisfy 

0<Xo<2^C\I\\\M\\l2, (13) 
0< A, < -^^(7-l)(72 + i07 + l)7-^ J = 1,2,3,... (14) 

To start proving these claims, let 

X,^\\M,\\L2^n,)^/\I^\ 

so that Mj/Xj satisfies the definition of an atom supported in the interval Ij. (This choice of Xj is 
positive except when Mj = 0, and in that exceptional case we may instead choose Xj to be any positive 
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number satisfying (13 1 or (14).) Thus the task is to estimate the norm of Mj and hence to estabhsh 
(131 and (14 1. First notice that for x G ilj, 



\M,ix)\ < \M{x)\ 



1 

W\ 



\M{y)\dy 



<|M(x)| + ^^||Af||i.(,,^.). 



(15) 
(16) 



Hence in particular, ||Mo||i2(Qg) < 2|| M||i2(Q^) by (16) and the triangle inequality. Thus Aq < 
2||A/||i2 y^(]T|, which is the desired estimate (13). 
Now consider j > I. Note that 



\Ij \ = and |rjj| = 2(7 - 1)7^-^7?. 



If X e rtj C (C/)'= then 



\M{x)\ < ^ , 



2RCi 



by the hypothesis ( |12[ ) with ?/q = 0, and so by ([15|), 

2RCi 



\M,{x)\ < 



Squaring and integrating shows that 



2RC4 , 

— r'^y 



(7 - i)7J--ii? cu 

2Cif R 1 



C (a;2 + 72j-ii?)- 



lAf, 



< 



(7-l)(7^ + 107 + l). 



Substituting this estimate into the definition of \j implies the bound (14 1. 
Step 2. From Step 1 we conclude that Y^Y=o belongs to , with norm 



4C4 /72 + 107 + 1 

~V 3(7-1) ■ 



j=o j=0 

Step 3. To construct the remaining parts of the decomposition of M, we let 

1 

V\ 



1, = 



n,- 



so that 



rrijlj = 



M{y) dy and 

M{y) dy)l,. 



3\ Jn, 



We will show X^j^o ^i^j belongs to . 



Let Tij = YlikLj "^fe- Then 



3=0 



by summation by parts and using that uq ~ J^M{y) = by hypothesis. Moreover, 



7-1 



7(7 + 1)^^--^^) 



(17) 
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is an atom supported in the interval Ij+i, because it has integral zero and because the disjointness of 
flj and flj+i implies that 

~ ~ /II \i/2 



^(7 + 1) 
7-1 



T 1-1/2 



Hence X^j^o '^i i belongs to H^, with atomic norm at most 



'7(7+1)^1 I 

' 3=0 



(18) 



Lastly, we observe that 



< 



ARC4 
4C4 



00 



R y 



dy by estimate (17) 



c 



-7 



Hence by (18 1 and a geometric sum, the atomic norm of X^j^o'^J i most 



AC4 rHn + i) 

C V (7-1)3 ■ 

Step 4. Clearly we may decompose M as 

00 00 

and so the bounds from Step 2 and Step 3 show that 



M\\hi < 2a 



+ — iy 3(7-1) 



'7^(7 + 1) 
(7-1)3 



The right side is minimal when 7 « 4.6. We choose 7 = 4, which lies close to the minimum point and 
gives a simple formula. The lemma follows immediately. □ 

Now we can deduce boundedness of Z on atoms. 

Lemma 3.4 (CZOs are bounded on atoms). Suppose Assumptions la, Ic, 2, 3 hold, and that a{x) is 
an atom. Let C > 3. Then Za G with 

\\Za\\Hi <2^\\Z\\l2^l-^+D{C)C3. 

Proof of Lemma 3.4 Since a is an atom, we have Za £ L^ f] L^ with J^{Za)(x) da; = by Assumption 
3, and so 

C,\I\ 



\{Za){x)\ < 



{x - yoY" 



for all X e (C-^)'^j by Lemma 3.2 (using Assumptions la, Ic and 2). Hence by Lemma 3.3 we conclude 
Za e iJi with 

4 



\Za\\Hi < 2^/C\r\\\Za\\L2 + —{V57 + ^320/3)0^. 



Now we finish the proof by substituting for C4 in terms of C3 (using the formula in Lemma 3.2) and 
recalling that ||a||i,2 < □ 

Next we develop a bound from L^ to weak-L^. The result is standard ^21i Theorem 8.2.1], and we 
include the lemma only to get a formula for the constant C5 . 
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Lemma 3.5 (CZOs map to weak-L^). Let ( > 3. If Assumptions la, Ic, 2 hold and f £ Ci , 
then 

>a}| < ^ll/IUi for all a > 0, (19) 

where C5 = •\/32(^ \\Z\\]^2_^]^2 + ^l^jCa. Hence Z extends to a linear operator from to weak-L^ that 
satisfies the same bound (19). 



Proof of Lemma 3. 5 First we show that if 6 G is supported in a bounded interval / centered at 
2/0 S K and h has integral 0, then 

2C 



\{Zh){x)\dx < 



Indeed, by ( 11 1 wc have 



when X ^ (^I. Integrating gives 



UZb)ix)\<C3 



UZb){x)\dx < C3 



C2-1 
\yo - y\ 

vo - y\ 



(20) 



\biy)\dy 



iJicir i^~y) 



2 dx\b{y)\dy. 



We evaluate the inner integral by changing variable like in the proof of Lemma 3.2 

HI- ^= f 

Y - yY 



{\w\>C} (w - Z)2 
2C|z| ^ 2C 



dw 



< 



since \z\ < 1 when y £ I. Hence we have proved (20 1. 

To prove the lemma, we now apply the proof of Theorem 8.2.1] to /, except changing the interval 
Q in that proof to J, and changing the interval Q* to (I, and using ( 20 ) as part of our estimate. One 
hence obtains the bound 



\{xeR:\{Zf)(x)\>a}\ < [8j\\Z\ 



+ 



c 

1^ e 



8C 



I/II 



where 7 > is a free parameter. (If we directly followed [21 then we would obtain 2C3 instead of C3, 
but we gain a factor of 2 by noting that the "bad" function in the Calderon-Zygmund decomposition 
has L-'^-norm at most 2||/||ii.) 

Minimizing over the choice of 7 yields the constant C5 stated in the lemma. □ 

Proof of Proposition [37l] Boundedness of Z on L^ for p € (1, 2), as stated in estimate ([t]), holds by 
combining the weak (1, 1) bound in Lemma |3.5| and the strong (2,2) bound in Assumption 2 with the 



explicit interpolation result in Proposition B.l (taking r — I there). 

Boundedness for p e (2, 00) as in ([s]) follows by duality, because Z* : L^ ^ L? satisfies Calderon- 
Zygmund Assumptions 1, 2 and 3 for the kernel K*{x,y) — K{y,x), except with the constants C2 and 
C3 interchanged. 

Now we prove Z is bounded on . Let / e Ll^ . Then / belongs to L^, and so Zf belongs to weak-L^ 
by Lemma 3.5 To obtain a formula for Z f , we consider an atomic decomposition / — X^j^i ^j'^j of 
/. Recall Zoj e with its norm bounded by an absolute constant, by Lemma 3.4 Hence the series 
^j^^j converges absolutely in . In fact, this series equals Zf a.e. by an argument that uses 
the weak-L^ result (see [5D1 p. 95]). Therefore the formula 

00 

Zf ^^XjZaj 
j=i 

holds, and holds for each choice of atomic decomposition. The norm estimate for Z in ([6]) now 
follows from the bound on atoms in Lemma 13.41 

Boundedness of Z on BMO as in ^ is immediate from duality and the boundedness of Z* on H^. 
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4. Wavelet frame operators are CZOs 

The wavelet frame operator ST is known to be a Calderon-Zygmund operator, under suitable hy- 
potheses on the synthesizer and analyzer. The existing proofs proceed in the spatial (or time) domain, 
using majorants: see |24l §5.6], [HI Chapter 7] or the explicit version in ^26, Lemma 3.1]. Un- 
fortunately, in practice such time domain estimates are too big (by an order of magnitude) to prove 
invertibility for examples such as the Mexican hat. We need better estimates that can capture cancela- 
tions in the analyzer and synthesizer. 

Theorem |4.6| in this section develops improved estimates on those Calderon-Zygmund constants of 
the wavelet frame operator, by working in the frequency domain and keeping ^ and (j) together there 
(rather than "splitting apart" ip and (fi in the time domain, as one does in the standard approach). Then 
later we use the theorem to prove our main result, in Section [5j 

4.1. Wavelet definitions. Recall the frame operator associated with synthesizer and analyzer 0: 

Formally, we may rewrite the frame operator as an integral operator 

{STf){x)^ f K{x,y)f{y)dy, 
Jr 

where the wavelet frame kernel is defined formally by 

K{x, y)=Bj2 {x)J~AF). X, 2/ e M. (21) 

Denote the portion of the kernel at level j = by 

Ko{x,y)^BY,Hx-Bk)(l){y-Bk). 

fcgZ 

4.2. Kernel kernel estimates in the spatial domain. We will prove absolute convergence of the 
series defining the kernel in the time domain, and of the series for its partial derivatives. Even though 
we ultimately want estimates in the frequency domain, the convergence of the kernel in the time domain 
ensures that the Calderon-Zygmund Assumptions 1 and 2 make sense. 

Lemma 4.1 (Convergence and differentiability of kernel at level 0). Assume ^ is hounded and \4>{x)\ < 
1/(1 -|-x^). Then the series 

Koix, y) = S ^ V(a; - Bk)^{y - Bk) 

is absolutely convergent, for each x,y € M.. Further: 

(i) If Ip and (j) are continuous, then so is Kq. 

(ii) If ip is continuously differentiable, ^' is bounded, and 4> is continuous, then the partial derivative 



= B^i)' {x - Bk)(j){y - Bk) 
feez 



exists and is continuous. 



The assumption \4'{x)\ < l/(l+a;^) has been chosen for its simplicity, and could certainly be weakened. 
It suffices to assume that (j) has a bounded, symmetric decreasing, integrable majorant, such as 1/(1 -|- 
jxl^"*""^). Similar comments apply to other lemmas below. 

Note. A formula for dKo/dy follows by interchanging the roles of and (/). 



Proof of Lemma \4-i\ For the absolute convergence, observe that 

\Ko{x,y)\ < (const.)||^||ic.^l/(l + (y-Bfc)2) 

fcgZ 

(i) Now suppose in addition that ip and (p are continuous. Then Kq{x, y) — > Kq{x, y) as [x, y) — )■ {x, y) 
by an application of dominated convergence. 
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(ii) By definition of tlie partial derivative 
dKo 



B lim ^ ^ [i}{x + h - Bk) ~ ij{x - Bk)] (f){y - Bk) 
fcez 

B^i^\x- Bk)<f){y - Bk) 



dx - - 

fcez 



using dominated convergence, since tlie difference quotient for "0 is bounded by HT/i'llioo. Continuity 
follows by part (i), since ip' is continuous. □ 

Now we show that the series for the full kernel K converges. 

Lemma 4.2 (Convergence and differentiability of kernel). Assume \'4>{x)\ and \(t>{x)\ are < 1/(1 + a;^). 
Then the series for K{x, y) converges absolutely whenever x ^ y, with 

\K(x,y)\ < B V \AYU(A^x - BkUiAJy - Bk)\ < (const.) , ^ , . 

^-^ \x — y\ 

Further: 

(i) Ifip and (j) are continuous, then so is K{x, y) wherever x ^ y. More precisely, the series K{x, y) = 
SjGZ l-^l'' 1^0 i-^'' X, A^y) converges uniformly on the set {{x,y) : \x ~ y\ > s}, for each e > 0, and hence 
K is continuous there. 

(ii) If ip is continuously difjerentiable, 4> is continuous, and \ip' (x)\ and \4>{x)\ are < 1/(1 + \x\^), 
then the partial derivative dK/dx exists by term-by-term differentiation wherever x ^ y, and dK/dx is 
continuous. This term-by-term derivative converges absolutely, and satisfies 

dK 

< H N I /I I 'J l-7/i' i AJ T — hih\rht AJii — Hk-W < I rnnsf I 

12 ■ 



dx 



<bY^ \A\^^^y{A^x- Bk)(f){A:iy - Bk)\ < {const.) j-^^^ 



Interchanging the roles of and (j) gives an analogous formula for dK/dy 



Proof of Lemma 4--S We adapt ideas from [24, Lemma 5.3.12 and §5.6]; or see |28| . 

Write r(x) = 1/(1 + x'^). For each k, the triangle inequality implies \x — y\ < \x — Bk\ + \y ~ Bk\, 
and so either — y| < — Bk\ or else i|x — y| < \y — Bk\. Hence we may "split apart" ip and 0, or 
rather the F-terms that bound them; that is, 



\Ko{x, y)| < C ^ F(x - Bk)r{y ~ Bk) 

< CT{^-\x ~ y\)Y,V{y - Bk) + i?F(^|x - y\) F(x - Bk) 
<CT{^\x-y\) 

where C is some constant (which changes from line to line, and depends on B). 
Summing the last estimate over j gives that 



\K{x,y)\<Y,\A\^T{A^x~y\/2) 
< (const.)-: 



F - y\ 

by Lemma |A.l| in the Appendix. 

(i) Let s > and |a; — y| > e. Since F is symmetric decreasing, the proof above gives 



\Kix,y)\ < J2 \Ay^{A^e/2) < (const.)- < oo. 



Hence the series K{x, y) — ^ \A\^ Kq{A^ x, A^y) converges uniformly on the set where |a; — y| > e. Since 
Kq is continuous by Lemma IXlJi), it follows that K is continuous when \x — y\ > e. 
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(ii) Write ri(a;) ~ l/(l + for the majorant of tp' and cj). By arguing as in the first part of this 
proof, one can show that the term-by-term derivative series B J^j kez — Bk)(t){A^y — Bk) 

converges absolutely on {x ^ ?/}, and is bounded by (const. )/|a; — yp. 

By Lemma |4.1[ ii) we know that this term- by-term derivative equals 

and that i^o is continuously differentiable with respect to x. This last series converges uniformly when 
|a; — ?/| > e, by arguments similar to part (i) above, and hence it converges to a continuous function. 
This continuous function is the derivative dK/dx. □ 

A function will possess majorants of the above type as soon as its Fourier transform possesses enough 
derivatives in L^, as the next lemma summarizes. 

Lemma 4.3 (Sufficient conditions for majorants). 

If ijj G W"^'^ then ip is continuous and \ipix)\ ^ + x'^). 
If ip G W^'^ then ip is continuous and \ip{x)\ < l/(l + ■ 

Proof of Lemma 4-3 From -0 G L"'^ we have that tp(x) = Jg^tp{^)e'^'^'^^^ is continuous (possibly after 
redefining ip on a set of measure zero). Integrating by parts twice gives that 

iP{x) ^ {2Trix)~^ f V^"(0e2"«"de 
Jr 

Hence 



|-2|1„/. "II 



^ 1 + x^ 



\^ix)\<mm[\mLi,\^^^rW 
li'lp £ W'^'^, then integrating by parts a third time yields a majorant that decays like \x\~^ at infinity. □ 



The hypotheses in the last lemma can be weakened to assume just that ip belongs to an inhomogeneous 
Besov space, Bf ^ or Bi ^ respectively. Those weaker assumptions do not appear useful for our work, 
though, because the Besov space norm seems too complicated for explicit estimations to be practical. 

4.3. Wavelet frame operator on — CZ Assumption 2. 

Proposition 4.4. Assume ip^cj) £ iP^ decay near the origin and infinity according to 

WOl<|,^,_v., |ei>i, (22) 

1^(01 < I, v., >i, (23) 

for some e > 0. Assume -0, </> G W^'^ . 

Then the frame operator ST satisfies C alder on- Zygmund Assumption 2: it is linear and bounded on 
, and if f £ has compact support and x ^ supp(/), then 

{STf){x)^ f K{x,y)f{y)dy. 
Jr 

Proof of Proposition \4-.4\ The frame operator ST is linear and bounded on L^, for example by [H Propo- 
sitions 6 and 7]. 

We know and (j) are < 1/(1-1- x^), by Lemma 4.3 and the hypothesis that 'ip,4> £ W'^'^. Thus by 



Lemma 4.2 the kernel satisfies \K{x, y)\ < l/\x — y\. If / € L has compact support and x ^ supp(/), 
then the function y i— l/\x — y\ is bounded on the support of /, and so by dominated convergence, 

K{x,y)f{y)dy^BY, f f{y)JjMy)dy^M^{STf){x). 



□ 
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4.4. Wavelet frame operator on weighted — CZ Assumption 3. 

Proposition 4.5. Assume </> S W^''^ with V'(O) — (/'(O) — 0, and that ijj' and (j)' decay near the origin 
and infinity according to ([T]). 

Then the frame operator ST satisfies Calderon-Zygmund Assumption 3: if f E has compact 
support and integral zero, then ST f € C\ with J^{STf){x) dx = 0, and (ST)* f G n with 
J^iSTrfdx^O. 



Proof of Proposition \4.^ If / G has compact support and integral zero, then its Fourier transform 
vanishes at the origin and belongs to the Sobolev space W^''^. Hence by boundedness of analysis and 
synthesis on the space of such "Littlewood-Paley functions" (as in [3J Propositions 5 and 6]), the wavelet 
series defining STf converges in the Littlewood-Paley space. Then the wavelet series converges also in 
L^, by an easy imbedding [21 formula (1)]. 

Integrating the series term-by-term gives that 



/ {STf){x)dx^ ^ B{f,(j)j^k)L^ / tpj^k{x)dx 



= 

because J^ip{x) dx — ^p{0) ~ 0. By interchanging the roles of (j) and ip we see the same conclusion holds 
for the L2_adjoint, since (ST)* = (S^T^)* = S^T^. □ 

4.5. Wavelet frame kernel — CZ Assumption 1. As explained at the beginning of the section, 
we need good, explicit estimates on the Calderon-Zygmund constants of the wavelet frame kernel. Our 
estimates involve the following frequency-domain norm quantities: 



ai(V,0)=^||V(-)'^(- + ZB-i)llLi 



a2{^P,cl,) = 2nJ2\\X{-m-)H- + lB-')\\L^ 
/ez 

a3(V,0) - 2nJ2Ui-)k- + IB-')W)\\l^ 



and 



ri(-0,0 


) 47r2 ^ 

/ez 


r2(V',0 


^ " 4^^ 

/ez 


T-3(V',0 


^ " 4^^ 



]"'IL. 



where X(^) = ^ denotes the identity function. 

Now we can estimate the Calderon-Zygmund constants Ci , C2 , C3 . 



Theorem 4.6 (CZ Assumption 1). Assume ^,(1) e L'^ with -0, 0, A-^, Ac?!) e W^^^ n W^'^ . 

Then the wavelet frame kernel K{x,y) in (21 1 satisfies Calderon-Zygmund Assumptions la, lb, Ic 
with the constants 



\A\ 

provided these constants are finite. 



^3(^, 0) = ^'^^y^l^^ ^a3(V'»r3(V,0)^ 
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(Aside. These constants blow up as \A\ — >■ 1. This apparent smgularity would disappear if we had 
multiplied the frame operator by log \A\, as discussed in the Introduction. Anyhow, the issue is irrelevant 
because we fix the dilation factor A in this paper.) 

The proof will build on off-diagonal decay lemmas for Kq and its derivatives. The first lemma explains 
the central idea: we invoke Poisson summation to transfer to the frequency domain, and then integrate 
by parts. 

Lemma 4.7 (Decay of Kq). Assume & L'^ with 'ip,(j> G W"^'^ C\ W"^'"^ , and that ai{ip,(p) < oo and 
Ti{ijj ^ <j)) < oo . Then for all x,y G M., 

n(V'» ■ 



|ifo(x,2/)| <min|ai(^,0),p^| 



X - y\ 

Proof of Lemma \4-. 7\ The hypotheses guarantee that tp and (j> are continuous and are < 1/(1 + a;^), by 
Lemma |4.3[ Hence the series for Ko{x,y) converges absolutely and yields a continuous function, by 
Lemma [4. II We have 

Koix, y)^BY, ^{Bk + x)0(Bfc + y) 

feez 

^Y,{^-^*^y^)i^B~^) (24) 

- e^'^*^^. By writing out the convo- 

(25) 



by the Poisson summation formula (justified below), where Ex{d 
lution, we find 



K^ix^y) ^Y.^-^^^'^'^'y I e2-«(--«)^(C)^(^ + IB-i)de 

This last series for Kq converges absolutely, with \KQ[x,y)\ < cri(V', 0) by definition of ai. 

Now we justify the use of Poisson summation. Let F{t) — Ko{x + t,y + t), so that F is continuous 
and B-periodic. The Fourier coefficients of F have £^-norm at most ai{ip,(t)) < oo, by the calculations 
above. Thus the Fourier series of is a continuous function, which therefore equals F{t) for every t. 
Choosing t = gives the Poisson summation formula as used in (24|, since F{0) = Ko{x, y). 
now to tl 

Ko{x,y) 



Returning now to the formula for KQ(x,y) in (25), integrating by parts twice gives that 

1 



y-g-2rfs-s f e2-»«^-2')[^(^)^(^ + ;B-i)]"dC, 



{2ni{x-y)y 

which implies that \Kq{x, y)\ < Ti{ip, <p)/\x — yp, and finishes the proof of the lemma. 



□ 



Lemma 4.8 (Decay of dKo/dx). Assume ?/','/'€ with e W^'^ andX^j;,(j} € W^'^ , and that 

a2{ip, 0) < oo and T2{ip, (j)) < oo. Then for all x,y £ M, 



dKo 



dx 



< min 



in|(T2(V', 



\x - y|3 . 



Proof of Lemma [ 



The hypotheses ^, Xip, (p S W'^'^ imply by Lemma 4.3 that ip, ip' and </> are contin- 
uous and are < 1/(1 + a;^). Thus the series for Kq and OKq/Ox converge absolutely and yield continuous 
functions, by Lemma |4.1[ 

Since OKq/Ox = B '^i^^j^ ip' {x—Bk)(j){y — Bk), the lemma follows by applying the proof of Lemma 4.7 
to tp' instead of ip, and integrating by parts three times instead of twice. 



Proof of Theorem 4.6, The hypotheses of the theorem ensure (by Lemma 4.3) that ip,ip' 
all < 1/(1 + |a;p). RecaU that the kernel 

i^(a;,y) = 5]|^pXo(A^x,^^y). 

has partial derivative 

\Ay A^ ^:^{A^x,A^y) 



are 



dx 



dx 
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by Lemma [42| ii). The first Calderon-Zygmmid estimate ^ follows by combining Lemmas 4.7 and A.l 



(in the Appendix), with z 



y and with Ci as defined in Theorem 4.6 Next 



dK 
dx 



< 



Co 



4 a;- 



V\ 



by combining Lemmas 4.8 and A.l with C2 as defined in Theorem 4.6 This derivative bound implies 



the second Calderon-Zygmund condition ([s]); for example, see jlJ, pp. 241-242]. The third estimate Q 
follows by swapping the roles of V' and cj). The theorem is now proved. 



5. Proof of Theorem wavelet invertibility on Hardy space, L'p and BMO 

Recall the function D{() defined for C > 1 in Section [2j and the constants Ci{ip,(f>) and so on 
introduced in Section 145] Define 



iVi(^, 0, C) = 2^\\ST\\l2^l^ + DiOC.iij, 0), 
N^iij, 0, C) = 2^\\ST\\l2^l- + D{0C2{i^, 0). 

Discarding the terms with D{Q gives an inequality on the norm: 

2\\ST\\l2^l2 <mm{N,{i^,^X),N^{ij,<f>,C)}- (26) 

Obviously, if iVi is finite for any value of C then it is finite for every value of and similarly for iVoo. 

Proposition 5.1 (Boundedness of the frame operator o n the Hardy space and BMO). Assume tp and 
(f) satisfy the regularity and decay hypotheses in Theorem \l.l\ that is, ip,(j),Xtp,X(f) G W^'^ O W^'^ and 
■!/'(0) = 0(0) = 0, along with estimate ([T]). Let > 3 and assume Ni{ip, <j), Q < 00 anrf A^oo(V'i '/'i C) < 00. 
Then the wavelet frame operator ST on extends to a bounded linear operator on the Hardy space 
and on BMO, with norm bounds: 

\\ST\\h^^h^ < A^i(^,0,C), 
\\ST\\bmo^bmo < Naoiip,(t), C)- 

Proof of Proposition \5. 1\ The frame operator ST satisfies Calderon-Zygmund Assumptions 1, 2 and 3, 
by Theorem |4.6[ Proposition |4.4| and Proposition |4.5[ respectively. Now Proposition |3.1| completes the 
proof. □ 



Defining the constants Mi and Mrx,- Next we define the constants in the statement of Theorem 1.1 
Assume tp and <j) satisfy the regularity and decay hypotheses of Theorem |1.1| (as in Proposition 5.1 ) 
Choose some other functions and that satisfy the same assumptions, and which in addition give 
perfect reconstruction on L^, meaning 

S'v,.?^. = Id : ^ (27) 

(Such functions certainly exist, for example, smooth orthonormal wavelets with sufficiently many van- 
ishing moments.) 
Define the constant 

Mp(V',(/)) = inf ^¥^(1/^-^*, 0,0 + infiVp(?/'*,0 -</-*,??), (28) 

C>3 77>3 

when p = 1,00. (This notation suppresses the dependence of Mp on our choices of and (j)^, but no 
harm will come from that omission.) 

Roughly speaking, if Mp is "small" then the perturbations tp — ip^, and 0—0* are "small" . Note that 
Theorem 1 1 . 1 1 assumes Mi and M^o are both rather small, in fact, less than 1 in magnitude. That means 
the perturbations are "small" , so that tp and are "close" to giving perfect reconstruction. 
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Proof of the main result — Theorem In view of the linearity of 5* with respect to the synthe- 
sizer -0, and linearity of T with respect to the perfect reconstruction in (271 implies a decomposition 

"ST — Id = S'^_^^T0 + S'^^T0_0_^ (29) 



valid on . Thus by Proposition 5.1 ST — Id extends to a bounded operator on with norm bound 

II 5r - id|| Hi^Hi < iVi - V'* , C) + A^i (0* , - 0* , v) 

for each Ci?? ^ 3. Taking the infimum of the right side with respect to C and 77 gives Mi{tp, and so 
we have 

||5T-Id|Ui^Hi <Mi(V'»<l, 
\\ST-ld\\L2^L- < ^Mi(0,0) < ^, 



by hypothesis and ( 26 ) . 

We conclude the frame operator ST is bijective on and on L^, with inverse given by the norm 
convergent Neumann series 

C30 

{ST)-^ = (Id - (Id - ST))-^ = ^(Id - STf. 

Notice this inverse operator is well defined on + L^, because if / G H then the inverse 
series X]fcLo(Id — ST)^f converges to the same limit function in as in L^, by considering pointwise 
convergence of a subsequence of partial sums. 

Since ST and {ST)~^ are bounded on and on L^, they extend to bounded operators on L^, 1 < 
p < 2, by interpolation (either the complex method [12, Corollary 1], or the real method 17, Theorem 
III. 6.1]). These extended operators remain inverses of each other on L^, by density of D L? . Thus 
ST is bijective on . 

Argue similarly for 2 < p < 00, using N^o and Mrx, and then interpolating between and BMO [T^ 
Corollary 2]. 



Note. An alternative to decomposition ( 29 ) is 

ST — Id = Sjf,Tcf,-ciy^ + S^-Tf,,Tcf,^ . 

If this formula is used then one obtains an alternative version of Theorem as one verifies by adapting 
the proof in the obvious manner. 



Comments on the Neumann series in L^. The quantities Mi and Moo determine the rate of 
convergence of the Neumann series for {ST)~^, in the proof above, for and BMO. To get convergence 
rates on one may interpolate Id — ST as follows. (This next material is not needed elsewhere in the 
paper.) Define 



Np{^,^X) = c{p){^\\ST\\l2^ 



(2/p)-l||_5,j.||2-(2/p) 



when 1 < p < 2 and let Mp{ip, (j)) be as in pS] ). Here c{p) is the constant from Section [sj Then 

\\ST\\l.^lp<Nj,[^,c^X) 

by Proposition 3.1 and so the decomposition (29) implies ||ST— Id||LP_j.Lp < Mp{tp, (f>). Hence if Mp < 1 
then the Neumann series for {ST)~^ converges with geometric rate Mp{ip,(l)). 

This approach breaks down for p near 1 because the interpolation constant c{p) blows up there, and 
so one gets Mp > 1. To obtain a convergence rate on the Neumann series for such p- values, one should 
instead interpolate powers of Id — ST, as explained in [3]. 
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Proof of Corollary |1.2[ Formally, one merely notes that bijectivity of ST on each function space in 
question (as given by Theorem implies surjectivity of 5, which means that each function / can be 
written as j, Cj^kipj.k for some coefficients Cj^k- 

To make this deduction rigorous, one needs S and T to be extended separately, whereas Theorem |1.1| 
says only that the composition ST extends from to be bounded on H^, BMO and L'p . Fortunately, 
the extension of analysis and synthesis separately to those function spaces (and indeed to a whole scale 
of Triebel-Lizorkin spaces) has been proved by Frazier and Jawerth Theorems 3.5, 3.7, and p. 81]. 
One can verify their hypotheses straightforwardly, given the assumptions in Theorem In particular 
they proved that T : — >■ /^'^ and S : /^'^ — ^ for a certain homogeneous Triebel-Lizorkin sequence 

space /p'^, and similarly for when p = 1 and BMO when p = oo. 
Thus for all / we have 

f = S{T{ST)-^f), 

and so / = with unconditional convergence, where the coefficient sequence {cj^k\ = 

T{{ST)-^f) belongs to the space J^-^. 

6. Example: Mexican hat wavelet expansions in Hardy space, LP and BMO 

The Mexican hat synthesizer ^Ij{x) = (1 — a;^)e~^'/^ is shown in Figure [l] along with its Fourier 
transform 

We assume dyadic dilations and unit translations, so that 

A = 2, B = l. 

We wish to apply our results to this example, in particular to get expansions of arbitrary f £ in 
terms of the Mexican hat system {ipj^k}- In order to apply our results, though, we must construct a 
suitable analyzer (p. 

1. First we will cut off the Mexican hat in the frequency domain, to obtain a band-limited synthesizer 
ip^,. Towards that end, define a "ramp" function 

{0 when ^ < 0, 

35^-* - 8Af + 70^^ - 20^^ when < ^ < 1, 
1 when ^ > 1, 

so that p G C3(M) and p'(^) = 140^^(1 - > for ^ G (0, 1). This formula for p' imphes that p{(,) 
increases from to 1 as ^ increases from to 1, and that 

PiO + /o(l - C) = 1 for all C G M. (30) 
Next define a cut-off function (see Figure [2]) by 

'p(6^-2) when^>0, 
when ^ < 0, 



<0 



so that K G C''(M) increases from to 1 as ^ increases from 1/3 to 1/2 and similarly as ^ decreases from 
-1/3 to -1/2. Further define 

V'* = (1 — k)^^. 

Obviously tpa, is band-limited, with -0* supported in the interval [—1/2,1/2]. Notice ip^. = tp on the 
interval [—1/3, 1/3], because k = there. 

2. Next we construct a (p such that analysis with (p followed by synthesis with ip^ gives perfect 
reconstruction on L^. The method depends on a dyadic partition of unity, constructed as follows. 
Start by defining a C'^-smooth "double bump" function supported in the interval [—1/3, 1/3]: 

rp(12^-l) whenO<C<i, 
m^\p{2-Q0 whenC>i, 
[Pi-O: whenC<0, 
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Figure 2. The cut-off function k and the double bump /3, in the frequency domain. 
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Figure 3. The Fourier transforms of the cut-off synthesizer, and of the analyzer. 



as plotted in Figure [2] We claim that /3 generates a dyadic partition of unity, with 

^ /3{2^0 = 1 for all C 7^ 0. 

Indeed, when < <^ < g we have (3(2^ S,) — for all j ^ 0,1 and the terms with j — 0,1 give 

m + mo = pm - 1) + p(2 - 120 = 1 



(31) 



by definition of [3 and ( 30 ). {Note. The double bump function /3 is a particular case of the "bell function" 
construction used for Meyer-Lemarie wavelets [24j.) 

We define the analyzer (j> by choosing its Fourier transform to be 

as shown in Figure [sj which is a kind of "band-limited reciprocal" of tp. Hence (p^p^ = (3 onR, because 
= on the interval [—1/3, 1/3], which contains the support of f3. Hence by the partition of unity 
property (311, we see -0* and p satisfy the "discrete Calderon condition" 



1, 



Since also and (p are supported in the interval [—1/2, 1/2] of length 1, the discrete Calderon condition 
implies that they give perfect reconstruction, meaning S^^T^ = Id on L^. (This fact, that band limitation 
together with discrete Calderon implies perfect reconstruction, can be found in [T71 Chapter 6], or else 
use [H formula (14)], for example.) 

Choose 0* — (j), so that ■0* and 0* give perfect reconstruction on L^. Then let 



Note that 'jl — nip. We plot /I(^) in Figure [sj 
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3. Next, for the synthesizer /i and analyzer (j) we estimate the quantities 
4>) = + 0M-)I1lM-i/34/3] < 0.000045 

a2(/x,(/>) = 2^^||/3(.)X(. + + 0M-)IUm-iAi/3] < 0.00022 

a3(M,0) = 2^ + 0M-)IIli[-iai/3] < 0.000067 

ti(m, 0) = ^ E II [/^(•)^(- + 0M-)]"|Li[_i/3,i/3] < 0-00086 

t2(m, ^ ^ ^ II + 0V^(- + 0M-)]"'ILi[_i/3,i/3] < 0-036 

-3(m, ^) = 4^ E II + 0M-)]"'|L.[_v3.i/3] < 0-014 

as we now justify. The first quantity cti was defined in Section [4. 5[ giving 

= eiim(-)0(-+oiili - Eii/^(- + O0(oiili. 

We may exclude / = from the sum because the supports of and cj) do not overlap (see Figure ^ . 
Further, when / 7^ we may replace 'p, (which equals K^p) with ip, because these functions agree outside 
the interval [—1/2,1/2]. Then we may substitute (j) — fSf-ip. These steps lead to the formula given 
above for (7i(/i, 0), which we have estimated numerically using Mathematica to obtain the upper bound 
0.000045. Wc proceed similarly for each of the other quantities (J2, 0-3, ti, r2, T3. 

From this numerical work we can estimate the Calderon-Zygmund constants in Theorem |4.6[ taking 
A = 2 there: 

Ciifi, (P) = 4VfTi(Ai,0)ri(/i,(/)) < 0.00079, 
40 , 

C2{^i. </>) = y < 0.088, 

40 



C3(m, = y ^a3(/i,(/.)r3(M,'/>)2 < 0.032. 

4. We have S'^^T^ = ST — Id on L^, since fi = tp — and S^^ T^ — Id on L'^. Hence 

\\S^T4l-^^l^ - \\ST-ld\\L2^L- < A 

where A denotes a Daubechies-type sufficient frame estimate; see [T, Theorem 1] with p — 2. The 
quantity A can be evaluated numerically, giving ||S'pT0||^2_^^2 < 0.00026. 

Using this fact and our estimates from Step 3, we calculate from the definitions of iVi and N^o that 

Mi{ij,cP) < iVi(^,0,9O) < 0.0075, 
Moo(V',0) < A^oo(/i>,180) < 0.011, 

where we chose the values of C = 90 and C = 180 to approximately minimize iVi and A^oo- (Recall also 
here that — </>* = 0, which eliminates the second term in the definitions of Mi and Moo-) Thus ST 

provides almost-perfect reconstruction on and BMO, with error of at most 1.1%. 

Hence the operator ST is bijective on and BMO, and on for each p e (1, 00), by Theorem |l.l[ 
We conclude from Corollary 1 1 . 2| that every function in can be expressed as a norm convergent series 
in terms of the Mexican hat wavelet system {ipj ^}. Thus our work provides a strong positive answer to 
Meyer's question in the Introduction. 
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Final comment on . Let us state some direct estimates for the frame operator on L^, in the Mexican 
hat example, even though we do not need them for our work: 

Mp{i}, (t>) < Npin, <j), 50) < 1 when 1.04 <p<2. 

When 1 < p < 1.04 we cannot argue directly that the Neumann series for (ST)^^ converges, because we 
cannot show Mp < 1 in that range. That explains why in the example above we use Theorem |1.1[ whose 
proof interpolates the invertibility property between and L^, rather than trying to prove directly 
that the Neumann series converges. 
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Appendix A. Geometric sums 

In order to extend our decay estimates from Kq to the full kernel K, earlier in the paper, we summed 
over dilation scales j S Z. Here we collect elementary estimates used in that task. 

Lemma A.l. Assume g -.M. ^ R and let u, r > 0. 

(i) If \g{z)\ < minjcr, r/|zp} for all z ^ 0, then 



(a) If \g{z)\ < min{(T, r/|z|'^} for all z ^ 0, then 



Proof of Lemma \A.l 



(i) Fix z ^ and let J be any integer. We split the sum into two pieces and estimate each piece 
using the assumption on g: 

,7-1 oo 

= j^^{\A\'-'<r + \A\-'\zrr) 

by the geometric series. Choose J to be the smallest integer satisfying jvll'^jzl > -^/r/cr, so that 

< -^/r/cr. (To motivate this choice, notice that a < T/|zp if and only if \z\ < y/rja.) 

Then 

^|^^g(A^z)| < j^{^\z\-'a+^\z\-\), 

from which part (i) of the lemma follows. 

(ii) Adapt part (i), except choose J to be the smallest integer satisfying |A|'^|z| > ^frja. □ 



Appendix B. Weak— strong interpolation between 1} and I? 

The Marcinkiewicz interpolation theorem implies in particular that if a sublinear operator satisfies 
weak-L"" and weak-L^ bounds, then strong-L^' bounds hold for p between r and 2. By strengthening 
the second hypothesis to strong-L^ and calling also on Riesz-Thorin interpolation (which requires the 
operator to be linear) , we will obtain an explicit estimate on the V norm such that equality holds when 
p=2. 
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ln(c(p)) 



Figure 4. The logarithm of the constant m Proposition B.l Note c(2) — 1 



The bound will involve the following constants, for 1 < r < p < 2: 

pr/fa+r) P{p + r){2-r) s(2-p)(p+r)/2p{p+r-pr) 

[p + r — pr)[p — r) / ' 

^ ^(2-p)(p+l)/2p 

P- 

Notice p = 2 gives c(2) — 1, and indeed c(2,r) = 1 for each r. At the other extreme of p-values, we see 
c{p, r) blows up like (p — r)^^/'' as p J, r, and c(p) blows up like {p — as p J, 1. See Figure [4] 

Proposition B.l. Let (X,^) and (Y,!/) be measure spaces, and 1 < r < 2. Assume Z is a linear 
operator from L'^{X) + L'^{X) to the space of measurable complex-valued functions on Y. 

If Z is weak type (r, r) and strong type (2, 2), then Z is strong type {p,p) for r < p < 2, with 



II 711 <- ^(r, 'r\ II 7||[(2/P)-l]/[(2/'-)-l] 



\z\ 



[(2/r)- 



(2/p)]/[(2/r)- 



(Equality holds for p = 2.) In particular, when r = 1 the conclusion says for 1 < p < 2 that 



\\Z\\lp(X)^Lp(Y) < c{p) \\Z\\^^/^x-^Xweak-L^(Y) \\^\\h{X)^L^Y)- 

A similar result holds if we assume strong type {q, q), but the proposition focuses on strong type (2, 2) 
because so many operators in harmonic analysis are bounded on L^. 

The case r = 1 of the proposition improves by a factor of about 8 (when p is close to 2) on the best 
bound we found in the literature, which is an exercise in the monograph by Grafakos [SHI Exercise 1.3.2]. 
We follow the same approach as Grafakos, except in the proof below we employ the harmonic mean of 
1 and p instead of the arithmetic mean; the harmonic mean yields simpler formulas. 



2-(2/p) 



Proof of Proposition \B. l\ Suppose r < q < p < 2. Then Marcinkiewicz interpolation applied with 
r < q < 2 (see [Ml Theorem 1.3.2]) implies boundness of Z on L"^, with 



|[(2/9)-l]/[(2/r)-l] 



|[(2/r)-(2/g)]/[(2/r) 



^q — r 2 — q 

Next, Riesz-Thorin interpolation applied with q < p < 2 (see Theorem 1.3.4]) says that 



(32) 



\Z\ 



Lp^Lp 



< \\Z\ 



[(2/p)-l]/[(2/g)-l] 



[(2/9)-(2/p)]/[(2/g)-l] 



We substitute the Marcinkiewicz bound ( 32 1 into this last Riesz-Thorin bound, giving that 

.[(2/p)-l]/(2-q) 



\z\ 



Lp^Lp 



< 
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1 - [r/q) {2/q) - 1 



(33) 



[(2/p)-l]/[(2/r)~l] 
I L^— >-wcak-L'" 



z 



[(2/r)-(2/p)]/[(2/r) 



We would like to choose q £ (r, p] to minimize the right side of bound ( 33 ) . That coefficient seems too 



complicated for analytical minimization to be feasible, but numerical work suggests that a good {i.e., 
within a factor of about 2 of being minimal) choice for q is the harmonic mean of r and p: 

1 lA 1. 
— -2^-r + -^- 



q 



p 



Substituting this choice of q into (33 1 yields the constant c{p, r) claimed in the proposition. 
Notice equality holds when p — 2, because equality holds in the Riesz-Thorin bound. 



□ 
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